arXiv: 1505.02517vl [physics.optics] 11 May 2015 


Impact of nonlinear loss on Stimulated Brillouin Scattering 

C. Wolff, 1 ’ 2 P. Gutsche, 2 M. J. Steel, 1 ’ 3 B. J. Eggleton, 1 ’ 4 and C. G. Poulton 1 ’ 2 

1 Centre for Ultrahigh bandwidth Devices for Optical Systems (CUDOS), 

2 School of Mathematical and Physical Sciences, University of Technology Sydney, NSW 2007, Australia 
3 MQ Photonics Research Centre, Department of Physics and Astronomy, 

Macquarie University Sydney, NSW 2109, Australia 
4 Institute of Photonics and Optical Science (IPOS), 

School of Physics, University of Sydney, NSW 2006, Australia 
(Dated: May 12, 2015) 

We study the impact of two-photon absorption (2PA) and fifth-order nonlinear loss such as 2PA- 
induced free-carrier absorption in semiconductors on the performance of Stimulated Brillouin Scat¬ 
tering devices. We formulate the equations of motion including effective loss coefficients, whose 
explicit expressions are provided for numerical evaluation in any waveguide geometry. We find that 
2PA results in a monotonic, algebraic relationship between amplification, waveguide length and 
pump power, whereas fifth-order losses lead to a non-monotonic relationship. We define a figure 
of merit for materials and waveguide designs in the presence of fifth-order losses. From this, we 
determine the optimal waveguide length for the case of 2PA alone and upper bounds for the total 
Stokes amplification for the case of 2PA as well as fifth-order losses. The analysis is performed 
analytically using a small-signal approximation and is compared to numerical solutions of the full 
nonlinear modal equations. 


I. INTRODUCTION 


In recent years the field of opto-mechanics has broad¬ 
ened from quantum-opto-mechanical research under¬ 
taken in high-Q resonators [l| to include the interaction 
of light with vibrations in high index-contrast optical 
waveguides [2]. The dominant opto-mechanical effect 
to occur in waveguiding geometries is Stimulated Bril¬ 
louin Scattering (SBS), which is the scattering of light 
from the travelling grating that is formed by an acoustic 
wave in the optical medium [3]. SBS was first proposed 
by Brillouin Q and subsequently observed in various sys¬ 
tems, ranging from first experiments with quartz |5], over 
the field of fibre optics, where it is well-known as a strong 
third-order nonlinearity || to more recent studies in on- 
ship waveguide such as chalcogenide rib waveguides 
This evolution has led naturally to the most recent inves¬ 
tigation of SBS in silicon nanowires UM based on ex¬ 
citing theoretical work [8| in the past few years. The very 
strong acousto-optical interaction that can be achieved in 
this system provides the potential to implement a num¬ 
ber of established SBS-applications in on-chip platforms; 
this includes novel light sources [12h-l3. non-reciprocal 
light propagation [15|, LL6| , slow light [17| , and signal pro¬ 
cessing in the context of microwave photonics |18l-[2Qj. 

Conventionally, the description of SBS is based on the 
approximation that SBS is by far the dominating nonlin¬ 
ear effect. A consequence of this approximation is that 
the total amplification of the Stokes wave along the wave¬ 
guide’s full length should be proportional to the power of 
the injected pump beam. The Stokes wave should then 
initially exhibit exponential growth until it starts to de¬ 
plete the pump; additional Stokes power can be gained 
by using a longer waveguide or by increasing the pump 
power. However, nonlinear loss is known to have rele¬ 
vant impact on the related effect of stimulated Raman- 


scattering 0 (especially in silicon photonics) and recent 
experiments in silicon waveguides at telecom wavelengths 
have also strongly suggested that nonlinear loss has an 
appreciable impact on the overall dynamic of the SBS 
process [ll|. In this context, third-order loss - i.e. two- 
photon absorption (2PA) - can be expected to have a 
qualitatively different effect on the SBS gain compared 
to fifth-order processes such as three-photon absorption 
(3PA) and 2PA-induced free carrier absorption (FCA). 
The latter is of considerable importance for semiconduc¬ 
tor waveguides, which are the best candidates for afford¬ 
able highly integrated optical circuits. Evaluating the 
impact of nonlinear loss terms on the SBS gain is critical 
for the design of future SBS-based devices. 

In this paper we study the impact of 2PA and 2PA- 
induced FCA on the performance of SBS. To this end, 
we derive nonlinear loss coefficients and solve the SBS- 
equations analytically within a small-signal approxima¬ 
tion. We thereby deviate from the Raman literature E3, 
which is necessary because SBS (in contrast to Ra¬ 
man scattering) is not always a forward-process, which 
would lead to directly integrable differential equations. 
Our analytical solutions provide strict upper bounds for 
key quantities such as output powers and amplification, 
which we express in terms of figures of merit for SBS in 
the presence of different types of nonlinear loss. For the 
case of negligible fifth-order losses (i.e. where FCA can 
be neglected), we find 

^ pA = F, (i) 

whereas the appropriate figure of merit for the case 
of dominant 2PA-induced FCA with a weak 2PA- 
perturbation is 

r-2/3 


‘Ffca = 


2v^7 ’ 


( 2 ) 
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where a, /?, 7 and T are the effective linear loss coef¬ 
ficient, effective 2PA-parameter, effective fifth-order loss 
(e.g. 2PA-induced FCA) parameter and the SBS-gain, re¬ 
spectively. These figures of merit must be greater than 1 
in order for the Stokes wave to be amplified. The applica¬ 
tion of these figures of merit can result in upper bounds 
for the Stokes amplification; an investigation of the spe¬ 
cific limits resulting from the presence of free carriers has 
been submitted separately as a rapid communication (23[. 

The layout of the manuscript is as follows: in Sec¬ 
tion an we state the preliminaries of our analysis, we 
introduce the relevant equations of motion, define and 
apply the small-signal approximation for the Stokes wave 
and state the expressions for the effective nonlinear loss 
coefficients for the case of intra-mode forward and back¬ 
ward SBS. In Section mu we solve the small-signal equa¬ 
tions analytically for systems that exhibit 2PA and lin¬ 
ear loss, for systems with 2PA-induced FCA with and 
without linear loss and we discuss the leading order per¬ 
turbative expression for case of linear loss and FCA with 
additional weak 2PA. Based on these analytical solutions 
we proceed in Section [IV| to derive figures of merit and 
present design guidelines for maximising SBS gain in ar¬ 
bitrary waveguide geometries. Conclusions and implica¬ 
tions are discussed in Section [VI Finally, we include three 
Appendices, in which we state the nonlinear coefficients 
for the more general case of inter-mode SBS as well as 


their derivations. 


II. PRELIMINARIES AND APPROXIMATIONS 

We consider the interaction between optical and acous¬ 
tic fields in waveguides with a material cross section that 
is invariant along the z -axis and supports both optical 
and acoustic guided modes. For the materials we as¬ 
sume the absence of magnetic response {fi r = 1), of ma¬ 
terial dispersion, of even-order nonlinearities (in particu¬ 
lar piezoelectricity) and of the Kerr effect. We explicitly 
include weak linear and odd-order nonlinear optical loss. 
The formulation is based on our previous work on the 
theory of SBS in integrated waveguides 0- We now in¬ 
troduce the equations and approximations that are solved 
throughout the remainder of this paper. We use SI-units 
throughout. 

A. Local acoustic approximation, power equations 

The starting point of this paper are the stationary cou¬ 
pled mode equations for the optical Stokes amplitude 
a^\z) at angular frequency cthe optical pump am¬ 
plitude a^ 2 \z) at angular frequency +f}, and 

the acoustic amplitude b(z): 


d z a (1) + (di + 4u|a (1) | 2 + 7m|a (1) | 4 ) a (1) = - (2/3 12 + 47n 2 |a (1) | 2 + 7i22|a (2) | 2 ) |a (2) | 2 a (1) 
d z a ^ + (d2 + ^22|ct.^ 2 ^| 2 + 7222|ct^| 4 ) a (2) =— (2^21+47221 |a (2) | 2 +7211 |a (1) | 2 ) | a ( 1 )| 2 a ( 2 ) 

d z b + a b b = - 

rb 


i uWQ 
p{ 1) 

ioj^Q 
Pi 2 ) 


a (2) &*, 

* 

-a™b. 


( 3 ) 

( 4 ) 

( 5 ) 


where the acousto-optic coupling parameter Q and the 
respective (signed) modal power fluxes V ( 1 ’ 2 ) appear as 
introduced in Ref. 0 - Note that the modal powers can 
be negative for modes that travel backwards (i.e. in the 
negative ^-direction) and that the acousto-optic coupling 
is real-valued in the absence of loss. The quantities cb, 
Pij and 7 ijk are the modal linear and nonlinear loss coef¬ 
ficients as derived in the Appendices of this manuscript. 
They are explicitly stated in Eqs. (|AldA 3 |) . Note the 
nontrivial factors 2 and 4 in front of the terms involving 
Pi 2, P21 1 7112 and 7221- 


As a first simplification, we assume that the acoustic 
decay length is much smaller than the length scale on 
which the optical envelopes vary. This is a valid assump¬ 
tion for SBS in fibres and in moderately long (mm-scale) 
integrated waveguides [7]. As a result, we can approxi¬ 


mate the acoustic profile: 


if }Q[a^ (; z)]*a^ (z) 
OLbPb 


(6) 


Under the assumption of phase matching and in conjunc¬ 
tion with the approximation a/ 1 ) « = cj, this allows 

us to eliminate the acoustic envelope in the equations of 
motion by inserting the substitutions 


vw>, 

with 

~ _ iAV\Q^_ 

vmp(2)v b a b 


( 7 ) 

(8) 

( 9 ) 
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into Eqs. da ni). The coupling parameter ceases to be along the waveguide: 
real-valued if phase-matching is broken I2ll . 

p(l,2) = ^(l,2)| 0 (l,2)|2 j (10) 

d z P {1 ’ 2) =2P (1 ’ 2) [a (1 ’ 2) ]* <9 z a (1 ’ 2) . (11) 


Next, we introduce the expression for the power fluxes 
in the Stokes and the pump modes and their derivatives 

I 


By inserting these expressions into Eqs. (l3l4l) . we obtain 
the equations of motion for the optical powers in a local 
acoustic approximation: 


d z P& + («i + p u P {1) + 7ni [P (1) ] 2 ) P {1) = - (2/3 12 - T + 47n 2 P (1) + 712 2 P m ) P (1) P (2) (12) 

d z pW + (a 2 + p 22 p( 2) + 7222 [p( 2 )] 2 ) p(2) = - ( 2^21 + r + 4 722 iP( 2 ) + 7 2iiP (1) ) P (1) p (2) , (13) 

with the power-related coefficients 

«.=2 /%=f§f. ijfjsj, r =2»{f). (14) 


B. Small-signal approximation 

We now introduce the central approximation of our 
work: We assume that at every position z inside the 
waveguide, the Stokes wave is much weaker than the 
pump wave: 

|P (1) (z)| < |P (2) (P)|. (15) 

As a consequence, a number of terms can be dropped 
from the Eqs. : all terms that are of at least sec¬ 

ond order in P W in Eq. and all terms that involve 
P W in Eq. (fl3|) . As a result, we obtain a simpler set 
of differential equations; one separable equation and one 
which depends on only a single quantity: 

d z P {1) = - ai P (1) 

+ (r - 2 / 3 12 - 7122 P( 2 ))P( 2 ) P (1 \ (16) 

d z P {2) = - (a 2 + /3 22 P (2) + 7222 [P (2) ] 2 )P (2) . (17) 

At first, this may seem to be a severe approximation. 
However, it can be motivated as follows: First, the ap¬ 
proximation leads to a distribution of the Stokes power 
along the waveguide that is strictly proportional to the 
injected Stokes power. This means that the SBS-active 
waveguide acts as a linear amplifier for the Stokes signal 
and only in this situation can the waveguide be expressed 
by an amplification factor that does not explicitly depend 
on the Stokes power. More importantly, this provides an 
upper bound for the amplification that is realisable for 
a certain set of parameters a,/?, 7 and T, because ev¬ 
ery term we neglect in the small-signal approximation 
introduces further loss. In other words, the Stokes am¬ 
plification predicted within this approximation is a up¬ 
per bound for the amplification that can be observed in 


reality with a finite Stokes input power. Stokes amplifi¬ 
cation can be observable in an experiment if and only if 
the small-signal approximation predicts it. 


C. Special cases: intra-mode SBS 

The expressions derived thus far apply to both inter¬ 
mode and intra-mode SBS, both in forward and back¬ 
ward configuration determined by the sign of the modal 
normalisation power V W of the Stokes wave (positive 
for forward scattering, negative for backward scattering). 
We now restrict ourselves to the more common case of 
interaction within the same branch of the optical disper¬ 
sion relation, i.e. we assume either intra-mode forward 
SBS (FSBS) or backward SBS (BSBS). The analysis pre¬ 
sented in the remainder of this paper can be carried out 
without this simplification and yields results of the same 
form, yet less transparent due to the required mode la¬ 
bels. In the case of FSBS, the Stokes and pump mode 
and respective normalisation power are identical: 

e (1) (x, y) =e (2) (x, y) = e(x, y) (18) 

-p(i) = -p(2) = p (1.9) 

If this is inserted in the expressions for the loss coeffi¬ 
cients derived in the Appendices, we find that all coeffi- 
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cients of equal order are identical: 


OL =OL\ = (12 

= -JT J d2r |e| 2 9{£r}, 

ft =Pll = @12 = @21 = @22 

J d 2 r ( 2 |e • e *| 2 + |e • e| 2 )S 2PA , 

7 =7111 : 7112 : 7122 = 7211 = 7221 = 7222 

2 


p 3 


/ 


d 2 r lei- 


|e • e | 2 + 2 (|e | 2 ) 2 


-iFCA 


( 20 ) 

( 21 ) 

( 22 ) 


where E 2PA and E FCA are nonlinear conductivities of the 
material associated with 2 PA and 2PA-induced FCA (see 
the Appendices for explicit expressions for values for the 
case of silicon). In the case of BSBS, the Stokes and 
pump mode are identical up to complex conjugation and 
the modal normalisation powers differ only in their sign: 


For the case of forward SBS with P( 1,F ) > 0 , only the 
sign of the right hand side of Eq. ([28]) flips: 

<9 2 P (1 ’ f) = - aP (1 ’ F) 

+ (r-2/3- 7 P (2) )P (2) P (1 ’ F) , (31) 


As a result, the forward SBS solution is exactly the in¬ 
verse of the corresponding backward SBS solution apart 
from a constant factor C: 


P^ F \z) 


C 

po^yyy- 


(32) 


The total forward amplification is thus identical to the 
backward Stokes amplification: 


•Tib (-0 — lOlogio 


-pPP)(L)' 

_p(i,F)( 0 )_ 


A%(L). 


(33) 


e ( 1 ) (x,y) =[e ( 2 ) ]*(x,y) = e(x,y), (23) 

pW = - p( 2 ) = _p. (24) 

As it turns out, the resulting loss coefficients like the 
modal powers are identical up to a sign. We find for the 
loss coefficients: 


oq = 

- a, 

0 L 2 

=a, 

(25) 

@H = 

-A 

@2 i 


(26) 

7 lij = 

-7, 

I2ij 

=7, 

(27) 


for all i and j, where the expressions for <a, @ and 7 are 
the same as given for the case of FSBS. 

III. ANALYTICAL SOLUTIONS 


For this reason, we can focus on either of the two types of 
SBS whenever we employ the small-signal approximation. 
This is the case for the remainder of this paper except 
for the discussion in Section mca where we compare 
our analytical expressions to numerical solutions of the 
large-signal Eqs. mm- We choose to base the following 
discussions on the equations for FSBS and will omit the 
label F where possible: 

<9 Z P (1) = - aP (1) + (T - 2/3 - 7 p( 2 ))p( 2 )p( 1 ) j ( 34 ) 
<9 z P ( 2) = - (a + /3P (2) + 7 [p( 2 )] 2 )p (2) . ( 35 ) 

The general approach is as follows: we first solve the 
equation for the pump field, then insert that solution into 
the equation for Stokes field and integrate the resulting 
differential equation. 


In this section, we present and discuss the analytical 
solutions to the small-signal equations for several experi¬ 
mentally relevant combinations of loss mechanisms. As a 
first step, we restate the equations of motion for the case 
of backward SBS. For convenience, we prefer the Stokes 
power to be positive definite, i.e. p(b B ) = P|a^^| 2 , lead¬ 
ing to the equations: 

d z P {1 ’ B) =aP (h b)_ 

(r —2/3 — 7 P ( 2 ’ B) )P ( 2 ) P (I ’ B) , (28) 

d z P {2) = - ( a + /3P (2) + 7 [P ( 2 ) ] 2 )P (2) . (29) 


The waveguide is assumed to extend from 2 = 0 to z = 
L. SBS-gain is usually specified in units of decibel per 
unit length of waveguide. However, this assumes that 
the Stokes wave grows exponentially, which is not the 
case in the presence of nonlinear loss. Thus, the natural 
quantity to compare is the total Stokes amplification in 
decibel over a waveguide of finite length L: 


*4!ib(^) — 10 l°gio 


-p(bB)( 0 )- 

_P(bB)(p)_ 


(30) 


A. SBS with 2PA and linear loss 

First, we investigate the case that the material does 
not exhibit fifth-order loss. This applies for example to 
glasses, including the chalcogenide rib waveguides that 
are the basis of recent on-chip SBS activities 0 , or 
to semiconductors illuminated below the half-bandgap. 
This case implies {<a,/d} 7 ^ 0; 7 = 0. The resulting equa¬ 
tion for the pump power is: 

a z P (2) + aP (2) + /3[P (2) ] 2 = 0. (36) 

This is an equation of the Riccati type and can therefore 
be solved using the ansatz P^ 2 \z) = y'/(()y), resulting 
in the solution 


P {2 \z) 
where ( 


<aexp(— az) 

/3[C - exp(-az)] ’ 

a 


=1 + 




( 2 ) 


0 


(37) 

(38) 
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normalized waveguide length (L / a) 


FIG. 1. Stokes amplification in decibel of a waveguide that is 
subject to linear loss and 2PA as a function of the pump power 
level and the waveguide length in natural units for a SBS gain 
ofT = 20 f3. Note that for every fixed pump power above a cer¬ 
tain threshold, the amplification grows with increasing wave¬ 
guide length, reaches a maximum and then drops off again. 
These maxima indicate optimal waveguide lengths and are 
connected with a solid black line. 


parametrises the input pump power. Next, we insert this 
result in the Stokes equation: 

P (1) , (39) 

which can be readily solved: 


d z P w = - ■ 


1 - 


(T — 2/3) exp(-az) 


(3[( - exp (-az)] 


d z P t 1 ) T — 2/3 a exp (—az) 

- = — - • -• 

Pt 1 ) [3 £ — exp (—az ) 5 

P^\z) =S 2PA exp(-az) [£ — exp(— az)]^ P ^~ 2 ^ , 


(40) 

(41) 


where 

S 2PA = p o (i) (c _ 1)( 2-r/« (42) 

parametrises the input Stokes power = P^(0). 

From Eq. m, we can readily derive the Stokes amplifi¬ 
cation in the presence of 2 PA and linear loss: 

A% A+Un (L) -aL + (T- 2/?) x (43) 



This result comprises two parts: First the exponential 
decay of the Stokes power due to the linear loss aL. Sec¬ 
ond, the effect of SBS and 2 PA. The logarithm is always 
non-zero and can be interpreted so that 2 PA and lin¬ 
ear loss reduce the SBS-effective waveguide length and 
SBS-effective pump power; it grows monotonically as a 


function of the waveguide length L and the pump power 
Pq , but saturates for aL 1 while being unbounded 
with respect to Pq 2 \ The prefactor (T — 2/3 ) finally is 
the central result of this analysis. It states that the SBS- 
gain of a waveguide with 2 PA present is effectively re¬ 
duced by twice the 2PA-coefRcient: T e q = T — 2/3. Thus, 
SBS can only be observed if T/(2(3) > 1. Therefore, 
although any level of total Stokes amplification can be 
obtained by increasing the pump power, increasing the 
waveguide length much beyond 1/a is not useful. To 
illustrate this, the general dependence of Al^ A+lm ^ for 
an arbitrary chosen ratio between SBS gain and 2 PA- 
coefficient of T = 20/3 is depicted in Fig. [1] over a wide 
range of waveguide lengths and pump powers. 

To conclude, we consider the case of vanishing linear 
loss a —>> 0. To this end, we approximate exp (—aL) « 
1 — aL in Eq. (ITT]) and Eq. (H3|) to find: 



with the corresponding Stokes amplitude 


P^\z) =P 0 (1) 


1 + z(3Pk 2) 


07/3—2) 


(45) 


We state this result to point out that in the case of 2 PA as 
the only loss mechanism, the solution is algebraic rather 
than transcendental and that any arbitrary exponent can 
be realised. For example, it is possible to connect the 
Stokes power linearly or quadratically to the pump power 
by choosing T = 3/3 or T = 4/3, respectively. We cannot 
think of any useful application for this analog computa¬ 
tion, but find it an amusing curiosity. 


B. SBS with FCA alone 


Next, we focus on the case {cq /?} = 0; 7 7 ^ 0, which de¬ 
scribes a system subject only to 2 PA-induced free carrier 
absorption, but neither linear loss nor 2 PA itself: 2 PA is 
assumed to generate carriers but not absorb a noticeable 
amount of energy itself. The following discussion is in¬ 
tended as a preparation for the discussion of 2 PA-induced 
FCA in combination with linear loss in Section IIII Cl and 
for the qualitative discussion in Section UlI D1 The results 
can also be directly applied to situations where both a 
and /? are very small, e.g. for very good semiconductor 
waveguides. 

Within this section, the pump power satisfies the equa¬ 
tion 


d z p^ 

which has the solution 

P {2 \z) 

where zq 


-i[p (2) ]\ 

1 

a /27 (z + zoY 

1 


(46) 

(47) 

(48) 
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FIG. 2. Solution to the problem of SBS with free carrier 
absorption in the forward SBS (top panel) and backward SBS 
(bottom panel) configuration for the solution normalisation 
constant set to *S FCA = VX. Note the respective extrema of 
the solutions at z + zq = A/2. 


again parametrises the input pump power. Note that 
negative values for zq are unphysical, because 7 is pos¬ 
itive. It is convenient to express lengths and powers in 
terms of the respective natural units of A (length) and II 
(power): 


a =7/r 2 , n =r/ 7 . ( 49 ) 

Next, we insert Eq. (1471) into Eq. (1341) to obtain: 


d z p « _ i r 

p (1 ) 2(z + 2 0 ) 7271 +^)’ 


(50) 



10 20 30 40 50 60 70 80 90 100 

normalized waveguide length (l_F/y) 


FIG. 3. Stokes amplification in decibel of a waveguide whose 
only loss mechanism is 2PA-induced free carrier absorption 
(2PA itself neglected) as a function of the pump power level 
and the waveguide length in natural units. Please note that 
for every fixed waveguide length, the amplification grows with 
increasing pump power, reaches a maximum and then drops 
off again. These maxima correspond to optimal pump power 
levels and are connected with a solid black line. 


waveguide and grows towards both ends. Furthermore, 
for very high pump levels, the Stokes amplitude at the 
output is lower than at the input, which means that any 
SBS-gain inside the waveguide is destroyed by free carrier 
absorption if the pump is too strong for a given waveguide 
length. 

The Stokes amplification of the waveguide can be de¬ 
rived from Eq. we find 


The solution to this equation is 

r\/2 (z + z 0 )/ 7 


cFCA 

P w (z) =—A=== exp 


y/Z -\- Z 0 

with the normalisation constant 

<5 FCA =P 0 (1) v^oexp -r72z 0 /7 . 


(51) 


(52) 


This solution Eq. (j5Tj) is plotted in the top panel of Fig. [2] 
for S FCA = VX directly above the inverse of Eq. 
i.e. the solution to the BSBS problem. The first feature 
of Eq. (I5T1) is that the general shape of the solution is 
universal for SBS-waveguides with free carrier absorp¬ 
tion. Within this plot the waveguide corresponds a win¬ 
dow starting at ^o/A (defined by the pump power) and 
of length LI A. Increasing the pump power simply moves 
the window to the left of the plot. Variations in the SBS- 
gain, the FCA-coefRcient and the injected Stokes power 
only rescale the plot axes. The second crucial feature 
of the solution is the extremum at z + zq = A/2. This 
means that in a strongly pumped waveguide the Stokes 
amplitude assumes its minimum somewhere inside the 


V / 27L(P 0 (2) )2 + 1-1 

-Iln[l + 27i(P 0 (2) ) 2 ]}- (53) 

This function is shown in Fig. [3] for a wide range of wave¬ 
guide lengths L and pump powers . As in Fig. [2j 
this plot has been made universal for any combination 
of SBS and FCA coefficients. As a consequence of the 
non-monotonic nature of Eq. m, the Stokes amplifica¬ 
tion A assumes a maximum at a specific pump power 
for every given waveguide length. These optimal pump 
powers can be computed by solving the equation 


a fca 

Aib 


(L) = 


io f r 

In 10 \ ry pA 


flo(2) U V 

L °- r 0 J L 

with L kept fixed. Since a closed analytical solution to 
this cannot be found, we use Newton’s method to find the 
zeros. However, it can be shown that the optimum pump 
powers always lies inside the interval 1/2 < 7 _1 TPq 2 ) < 
1. The numerically determined optimal pump powers are 
highlighted in Fig. [3] with a solid black line. 
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C. SBS with FCA and linear loss 

After the discussion of FCA as the only loss mecha¬ 
nism, we now add linear loss, i.e. we study the situation 
/3 = 0; {<a, 7 } 7 ^ 0. This is the most general case in¬ 
volving FCA that still can be solved analytically. In Sec¬ 
tion mna we present a perturbative treatment of weak 
2 PA alongside strong linear loss and FCA based on the 
expressions derived in this section. 

As before, we start with solving the equation for the 
pump power along the waveguide: 

d z P (2) = -aP (2) - 7 [P^] 3 . (55) 

This can be transformed into a Riccati equation via a 
substitution of the type u = [P^] 2 . The closed solution 
is: 


p (2) 0) = 


p 0 (2) v^ 


( 7 [P 0 (2) ] 2 + a) exp (2 az) - 7 [P 0 (2) ] 2 


(56) 


Next, we solve for the Stokes power along the waveguide 
by integrating the equation 


d z P {1) 

p(i) 


= -a + TP^ -7[p( 2 )] 2 , 


to obtain 


p(*)(z) =S FGA+lin p( 2 \z)x 


exp 



\fa 

_V7P (2) (z)_ 


(57) 


(58) 


where the input Stokes power enters via the normalisa¬ 
tion constant 


gFCA+lin 


>( 1 ) 


>( 2 ) 


exp 



\fa 


• (59) 


The spatial evolution of the Stokes wave is qualitatively 
similar to the one depicted in Fig. O By taking the 
decadic logarithm of Eq. (EH) including Eq. , one 
can directly obtain the explicit relationship between the 
total Stokes amplification ^4dB A+lm an d waveguide 
length and the pump power. This expression is long and 
convoluted and does not provide any additional insight, 
so we forgo showing it here. As before, it is very useful to 
introduce the problem-specific units of length and power 
stated in Eq. (03) as the natural unit system. Conse¬ 
quently, the linear loss is best expressed in natural units 
of inverse length: 


T = cryr -2 . (60) 

If this normalised linear loss vanishes, the plot of 
AdB A+lm identical to the previous result Fig. [3j As 
the linear loss increases, the total amplification decays 
not only for increasing pump power but also for increas¬ 
ing waveguide length, leading to a well defined maximal 
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FIG. 4. Stokes amplification in decibel of a waveguide that is 
subject to linear loss and 2PA-induced free carrier absorption 
(2PA itself neglected) as a function of the pump power level 
and the waveguide length in natural units for a normalised 
linear loss T = ajF~ 2 =0.03. To illustrate the impact of 
linear loss, we adopted the colorbar from Fig. [3] There is one 
optimal combination of waveguide length L^ opt ^ « GOyr -2 
and pump power p(°pt) ~ p/y that leads to the maximally 
realisable Stokes amplification of about 15 dB for this value 
of T, which corresponds to a figure of merit (see Section HVBI 
for details) of Pfca ~ 2.9. 


obtainable amplification for each value T. The maximal 
amplification is obtained for exactly one pair of optimal 
pump power p(°P t ) and optimal waveguide length l/ opt ). 
This is illustrated in Fig. 0] for the case of T = 0.03. We 
will resume this topic in Section IIVB1 


D. SBS with FCA, linear loss and weak 2PA 

Finally, we turn towards the case that all three loss 
mechanisms (linear loss, 2 PA and FCA) are present. In 
this case, the equation for P^ 2 \z) can only be derived 
implicitly, i.e. in the form z(P^), but not explicitly. Un¬ 
fortunately, this makes it impossible to solve the general 
case exactly. One therefore has to resort to approximate 
solutions derived from numerical calculations or pertur¬ 
bation theory in the parameter /?. To zeroth order in /?, 
the pump profile P^ (z) can be adopted from the 2 PA- 
free case Eq. (j47j). This is because we assume that the 
impact of 2 PA on the pump is of the same order as the 
(neglected) pump-depletion due to SBS. Under this as¬ 
sumption, the SBS-gain T is simply reduced by twice the 
2 PA-coefficient: 


r eff = r - 2 / 3 . (6i) 

All solutions from Section IIII Cl apply directly with this 
substitution. 

The next perturbation order would involve a first or¬ 
der correction to the pump profile P^ (z). While this can 
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be found, the subsequent equation for the Stokes ampli¬ 
tude P^\z) and as a consequence the equation for the 
total Stokes amplification involve an integral that does 
not have a closed solution. For this reason, we adhere to 
the zeroth order approximation and refer the reader to 
numerical solutions whenever more accurate results are 
required. 


This means that the optimal waveguide length depends 
logarithmically on the pump power, but will always re¬ 
main in the vicinity of aT 1 . The power-dependence of 
the optimal waveguide length for the case of J~ 2 pa = 10 
is shown in Fig. [l] as a black solid line. 

B. SBS with FCA, linear loss and weak 2PA 


IV. DESIGN GUIDELINES AND FIGURES OF 
MERIT 

In this section, we provide guidelines for the optimal 
choice of waveguide length and pump power of waveg¬ 
uides whose SBS gain and loss coefficients have either 
been measured or computed numerically using the ex¬ 
pressions in Ref. sa and the Appendices of this paper. 
We furthermore provide figures of merit that express the 
suitability of any given material or waveguide design for 
the purpose of SBS in a single number. 


A. SBS with 2PA and linear loss 


In the absence of fifth-order loss (especially 2PA- 
induced FCA), we can apply the results from Sec¬ 
tion IIII A1 This covers most centrosymmetric and amor¬ 
phous insulators including glasses. The main feature of 
Eq. (I43j) is that it is monotonic and unbounded with re¬ 
spect to the pump power. Thus, the SBS gain of any 
given waveguide can in principle be increased indefinitely 
by injecting a sufficiently strong pump, provided the SBS 
gain is sufficient to overcome the loss due to 2PA. In other 
words, any material or waveguide design is capable of ex¬ 
hibiting SBS if the figure of merit 

-E2 pa = prz (62) 


is greater than one. In every material or waveguide design 
with J- 2 pa < 1, SBS is quenched by 2PA. For T 2 PA = 1, 
SBS and 2PA cancel each other to leading order; higher 
order corrections predict a weak power-dependent decay 
of the Stokes wave in this regime. 

The total Stokes amplification does not grow indefi¬ 
nitely with respect to waveguide length, because the lin¬ 
ear term —olL in Eq. m at some point overcomes the 
saturating logarithm of the second term. The optimal 
waveguide length L( opt 0 [ s given by the condition 


&4^ A+liE 


dL 


= 0 , 


which can be evaluated exactly: 


£(° pt > =- In 
a 


p o 2 \ r -P) 

(/3P 0 (2) + «) 


Pq (2) >q// 3^ ln(2J r 2PA - 1) 
a 


(63) 


(64) 

(65) 


In most indirect semiconductors such as silicon and 
germanium, the long carrier lifetime leads to the ef¬ 
fect that the free carriers created by 2PA have a much 
stronger impact on the absorption of a quasi-CW light 
wave than the two-photon absorption itself. Although 
this effect can be reduced by extracting free carriers via 
an externally applied electric field, the loss due to 2PA- 
induced FCA will still surpass the 2PA itself in many 
situations and the results from Section IIII Dl can be ap¬ 
plied. 

In contrast to the previous section, the total Stokes 
amplification of a waveguide that experiences FCA is 
non-monotonic with power and is in fact bounded with 
respect to both the injected pump power and its length. 
For each set of SBS gain and loss parameters, there is 
exactly one choice of length and pump power that leads 
to the maximal gain, as can be seen in Fig. 01 

In any given material or waveguide design, a total 
Stokes amplification can be obtained only if the FCA 
figure of merit 

•Ffca = (66) 

2 v / «7 

is greater than 1. As shown in [23|, this quantity emerges 
naturally as a figure of merit from Eqs. cnn. Unlike the 
situation without FCA (see Section HV A[) , Pfca merit 
also implies an absolute maximum for the Stokes am¬ 
plification that can be obtained with a given waveguide 
design. It is attained for the optimal choice of wave¬ 
guide length and pump power, both of which depend on 
*Ffca- This interdependence is discussed in a separate 
paper [23]. 

In Fig. [5j we compare the maximally realisable Stokes 
amplification predicted by our small-signal approxima¬ 
tion with results obtained by numerically solving the 
large-signal equations Eqs. (mini) for the analytically 
predicted operating conditions and injected Stokes pow¬ 
ers between —60 dB and —20 dB of the natural power unit 
II both in forward and backward configuration. In the 
case of FSBS, we employed a simple 4-th order Runge- 
Kutta integrator. In the case of BSBS, we used the shoot¬ 
ing method based on a variable order Runge-Kutta that 
is optimised for stiff differential equations. It is clearly 
visible how the small-signal solution provides an upper 
bound for the realisable Stokes amplification. The differ¬ 
ence can be attributed to energy transfer from the pump 
to the Stokes (known as pump depletion in SBS with¬ 
out loss) and additional carrier generation due to the 
relatively high Stokes intensity. Finally, we note that 
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FIG. 5. Stokes amplification plots for linear loss, FCA and 
weak 2PA: Comparison of the maximally realisable ampli¬ 
fication within the small-signal approximation (black dash- 
dotted line) and corresponding maximal amplifications in 
FSBS (solid lines) and BSBS (dashed lines) configuration for 
finite Stokes input powers between — 60 dB and — 20 dB of the 
natural power unit II = r e ff/7 = (r — 2/3)/y. 

the optimal operating conditions (i.e. waveguide length 
and pump power) shift as the input Stokes level is in¬ 
creased and that this affects forward SBS and backward 
SBS quite differently. This is because the small-signal ap¬ 
proximation no longer applies over the full range of Tfca 
and input Stokes powers. In the case of BSBS, the am¬ 
plification maximum as a function of waveguide length 
and pump power remains very flat, which means that 
the analytical prediction is very viable. However, in the 
case of FSBS the amplification maximum becomes fairly 
sharp and the waveguide has to be shortened to obtain 
performance that is comparable to the BSBS amplifica¬ 
tion shown in Fig. [5j In other words: A forward SBS 
amplifier has a smaller dynamic range than an amplifier 
based on backward SBS. 


V. CONCLUSIONS 

In this paper, we have discussed the impact of non¬ 
linear optical loss on the process of SBS within a small- 


signal approximation. Based on analytical solutions to 
this problem, we have derived figures of merit that de¬ 
scribe the suitability of a material or a waveguide design 
for the case that fifth-order nonlinear loss effects such 
as 2PA-induced FCA can be neglected and the case that 
they cannot. In the former case, we find that, although 
third-order loss reduces the total Stokes wave amplifica¬ 
tion along the waveguide, the amplification is not funda¬ 
mentally restricted and can be increased indefinitely by 
increasing pump power until the weakest component in 
the optical circuit is permanently damaged. In the lat¬ 
ter case, however, the Stokes amplification due to SBS 
is overcome by FCA once some optimal pump power is 
exceeded. Thus, the total amplification is bounded to 
a value that is ultimately given by the figure of merit 
T fca = (T — 2(3)/^/4xry. This opens basically three 
possible routes to effective SBS-circuits in semiconduc¬ 
tor (especially silicon) photonics platforms. First, the 
upper amplification bound can be increased by reducing 
the linear loss a. This often is a challenging task. Sec¬ 
ond, the impact of FCA could be drastically reduced by 
removing the free carriers by means of an externally ap¬ 
plied electric field. The resulting increase in linear loss 
would be tolerable as long as the product is reduced. 
The third, and on a short time scale potentially most 
viable, solution would be to design future silicon pho¬ 
tonic circuits for SBS to operate at wavelengths below 
the 2PA-threshold, e.g. around 2400 nm for the case of 
silicon. This would eliminate not only 2PA (which itself 
is not problematic), but in particular the induced FCA 
and lead to higher acoustic quality factors due to the 
reduced acoustic frequency. 
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Appendix A: Expressions for loss coefficients 


In the following Appendices, we derive the loss coefficients that appear in Eq. (|3]) and Eq. (|4]). 
expressions are: 


The explicit 


~ _£o^_ 
* ~~ p(i) 

Pv ~ -p(i) 

_ _ 1 



(Al) 

(A2) 

(A3) 


where the indices i, j and k label the respective optical eigenmodes modes and can take the values 1 and 2. The 
Symbols E 2PA and E FCA are nonlinear conductivities associated with 2PA and 2PA-induced FCA and are expressed 
in more conventional terms in Eq. m and Eq. (El, respectively. The values 


E 2PA =5.5 • 10“ 16 Wm 3 V“ 4 , 


E FCA =6.0 • 10“ 28 Wm 5 V“ 6 , (A4) 


correspond to a bulk 2PA-coefficient of 5 • 10 -12 m/W and an electron scattering cross section of 1.45 • 10 -21 m 2 
together with a carrier life time of 10 ns. Those are typical literature values (24| for silicon at a vacuum wavelength 
of 1550 nm. 


Appendix B: Derivation of 2PA-terms 

In time-domain it is sometimes advantageous to represent lossy (i.e. dispersive) optical nonlinearities is via a non¬ 
linear current distribution [25j. The corresponding nonlinear conductivity to represent 2PA is related to the imaginary 
part of the third-order susceptibility: 


^ 2 PA 


= we 0 S{x (3) }. 


(Bl) 


Here, we neglected the tensorial nature of ;\' (3) in order to to improve the readability of the integrals (IATM31) . The 
generalisation to tensorial nonlinearities is a matter of book-keeping. With this, we find the real-valued time-domain 
current density J 2PA (t) = — E 2PA |E(t)| 2 E(t), where for the sake of brevity, we wrote v 2 = v • v and |v| 2 = v* • v 
for any vectorial quantity v. The electric field is formed by the interference between two optical eigenmodes: E (t) = 
a ( 1 ) e ( 1 ) + a ( 2 ) e ( 2 ) + c.c.. Within a coupled-mode theory, we need the projection of J 2PA on the optical eigenmodes, 


e.g. e 


(i). 


(e«|J 2PA )=-( J d 3 r [(e( 1 ))*-(a( 1 )e( 1 )+a( 2 )e( 2 )+c.c.)]E 2 PA |a( 1 )e( 1 )+a( 2 )e( 2 )+c.c.| 2 

= -6a« J d 3 r |eW| 2 E 2PA 


d 1 ) I 2 Ifd 1 ) I 


+ \CL 


( 2 ) | 2 |( 2 ) | 2 ‘ 


(flW)* J d 3 r ( e W .e( 1 ))*£2 p A (a (i))2 e (i) . e (D 
■2a< 2 > j d 3 r (e^)* ■ e^ 2PA a^ • (e( 2 ))* 

2 (a^)* J d 3 r (e« . e ( 2 ))*S 


2PA a^a^e^ ■ e< 2 > 


- oscillating terms ) 

' T&( 


=a^\a ^\ 2 


2a^\aW 


J d 2 r E 2PA (2|5( 1 ) • (S«)*| 2 + |S« -e^\ 2 ) 

j d 2 r E 2PA (|e^^ -e^ 2 + |e< x > • (e( 2 ))*| 2 + 


di)| 2 |A(2)i 


! ). 


(B2) 


(B3) 


(B4) 


The projection on the other mode e^ 2 ) follows from this by interchanging the mode index superscripts. 
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Appendix C: Derivation of FCA-terms 


As in Appendix[B]we intend to describe the loss via a time domain current density J FCA . It is due to the conductivity 
of a dilute plasma with carrier density 7V C , which can be described [26j using a Drude-Sommerfeld model with effective 
carrier mass m* and damping parameter 


J FCA (t) 


(u 2 3 + Wj)m* 


E(i). 


(Cl) 


Here, q is the elementary charge and the parameters ra*, uj a are material specific and therefore may depend on position 
within the waveguide. The carriers are created by two-photon absorption and destroyed via recombination on a time 
scale r c > 10ns. The order of this time constant is important because it covers around 100 acoustic cycles for typical 
Stokes shifts of 3 — 30GHz. Consequently, the carrier density is to be derived from the 2PA power loss averaged over 
a time-scale greater than one acoustic period: 


p2PA 

-Moss 


( E (t) • J 2PA (£) ) = £ 2PA ( |E(i)| 2 |E(i)| s 


= 2 £ 2 PA || oW | 4 [| eW - e ( 1 )| 2 + 2 (| e ( 1 )| 2 ) 2 


+ a 1 


(2) |4 


3 ( 2 ) . e ( 2 )| 2 + 2(| e ( 2 )| 2 ) 2 


+ 4 | a «| V 2 )| 2 [| e « • e ( 2 )| 2 + | e «| 2 | e ( 2 )| 2 + | e « • ( e < 2 >) 


* 12 


(C2) 


(C3) 


The carrier generation rate is this optical power loss divided by twice the photon energy, whereas the carrier recombi¬ 
nation rate is determined by N c itself and the carrier life time r c : dtN c = P } 2FA / (2huj) — Note that the carrier life 
typically is reduced near material interfaces and therefore position dependent. In our crude model we neglected carrier 
diffusion. If the field amplitudes vary on time scales larger than r c , the loss due to FCA is given by the equilibrium 
carrier density = r c P ^^/. This allows us to express the FCA-related nonlinear conductivity: 


. , „2v2PA 
^FCA _ WdT c q L 


2huj(uj 2 + ’ 


(C4) 


Within the context of our coupled mode theory, we require the projection of the FCA-current on the optical eigen- 
modes, e.g. the mode e^: 


(e (D|jFCA )= __L_^ J d 3 r [e (l)]* . [, 

aWlaWl 4 


a ( 1 ) e ( 1 ) + + c.c. 


UdQ T c 


_p2PA 

2^*^1088 


(C5) 


hbJ 
huj 

4a( 1 )|a( 1 )| 2 |a ( 2 )| 2 


2huj((j 2 + cj|)m 

J d 2 r .5«| 2 + 2(|g( 1 )| 2 ) 2 ' 

J d 2 r • 5^| 2 + 2(|g( 2 )| 2 ) 2 ' 

j d 2 r |e^)| 2 £ FCA jje^ • e^| 2 + |e^)| 2 |e^ 2 )| 2 + |e^ • (e^)*| 2 J. (C6) 
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